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ABSTRACT 

We set sensitive upper limits to the X-ray emission of four Type la supernovae (SNe la) using the 
Chandra X-ray Observatory. SN 2002bo, a normal, although reddened, nearby SN la, was observed 
9.3 days after explosion. For an absorbed, high temperature bremsstrahlung model the flux limits 
are 3.2 x 10"^^ ergs cni^^ s~^ (0.5-2 keV band) and 4.1 x 10~^^ ergs cm~^ s~^ (2-10 keV band). 
Using conservative model assumptions and a 10 km s~^ wind speed, we derive a mass loss rate of 
Af ~ 2 X 10~^ Mq yr~^, which is comparable to limits set by the non-detection of Ha lines from 
other SNe la. Two other objects, SN 2002ic and SN 2005gj , observed 260 and 80 days after explosion, 
respectively, are the only SNe la showing evidence for circumstellar interaction. The SN 2002ic X- 
ray flux upper limits are ~4 times below predictions of the interaction model currently favored to 
explain the bright optical emission. To resolve this discrepancy we invoke the mixing of cool dense 
ejecta fragments into the forward shock region, which produces increased X-ray absorption. A modest 
amount of mixing allows us to accommodate the Chandra upper limit. SN 2005gj is less well studied 
at this time. Assuming the same circumstellar environment as for SN 20021, the X-ray flux upper 
limits for SN 2005gj are ~4 times below the predictions, suggesting that mixing of cool ejecta into 
the forward shock has also occurred here. Our reanalysis of Swift and Chandra data on SN 2005ke 
does not confirm a previously reported X-ray detection. The host galaxies NGC 3190 (SN 2002bo) 
and NGC 1371 (SN 2005ke) each harbor a low luminosity (Lx ~ 3 — 4 x 10'*° ergs s^^) active nucleus 
in addition to wide-spread diffuse soft X-ray emission. 

Subject headings: galaxies: individual (NGC 3190, NGC 1371) — supernovae: general — supernovae: 
individual (SN 2002bo, SN 2002ic, SN 2005gj, SN 2005ke) 



1. INTRODUCTION 

Type la supernovae (SN la) are an important sub- 
class of supernova (SN) that are thought to arise from 
explosions of white dwarfs in binary systems, although 
the exact nature of their progenitor systems is largely 
a mystery. This ignorance is not due to a lack of ef- 
fort since SN la are the subject of intense scrutiny, not 
least of all because of their importance as cosmological 
probes. SN la have been u s ed to measure th e Hubble 
constant (|Hamuv et al.lll995HRiess et al.lll996f ) and have 
even give n strong evidence for a non-zero c osmological 
constant (|Riess et al.lll998l : iPerlmutter et al.lfr999,) . Our 
failure to identify SN la progenitors highlights a ma- 
jor gap in our understanding of stellar evolution in bi- 
nary systems, and presents a stumbling-block to under- 
standing the chemical evolution of galaxies, since la's are 
known to be efhcient producers of iron. In addition, un- 
til SN la progenitor systems are clearly identified, it will 
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remain difficult to convincingly eliminate evolutionary 
sources of systematic error in the observed trend of SN 
brightness with redshift. 

The current model framework for Type la SNe involves 
carbon deflagration/detonation in a white dwarf driven 
close to the Chandrasekhar limit by accretion. The most 
likely type of progenitor system (jBranch et al.|[l995l ) is a 
C-0 white dwarf accreting H/He-rich gas from a compan- 
ion, either from its wind or through Roche lobe overflow. 
Double degenerate scenarios with coalescing pairs of C-0 
white dwarfs are also possible, while sub-Chandras ekhar- 
mass white dwarfs are less likely (jBranch et al.| [l995). 
In the non-coalescing scenarios, there will be circum- 
stellar gas whose composition and geometry depend on 
the nature of the progenitor system. If the circumstel- 
lar medium (CSM) emits radiation, or absorbs radiation 
from the SN, this can be used to distinguish between 
types of progenitor systems. 

Thermal X-ray emission is expected to arise from the 
hot gas in the interaction region between the rapidly 
moving supernova ejecta and the wind from presupernova 
mass loss. Detecting this X-ray emission offers a direct 
and potentially quite sensitive probe of the amount of 
CSM. Because of its high sensitivity, broad bandwidth, 
unprecedented imaging capability, and ability to respond 
rapidly to a target of opportunity (ToO) request, the 
Chandra X-ray Observatory is the premier instrument 
to search for faint X-ray emission from transient point 
sources. 

In the following we present the analysis and interpre- 
tation of Chandra observations of four recent SNe la. 
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Two are fairly normal SNe la, while the others are pe- 
culiar cases showing strong Ha emission from circum- 
stellar interaction. The following section describes the 
observed targets, §3 presents the observations and tech- 
niques, §4 uses the upper limits on X-ray flux to constrain 
the nature of the CSM, and the final section concludes. 
In an Appendix we describe the X-ray properties of the 
host galaxies for the two nearby systems, NGC 3190 and 
NGC 1371. 

2. DESCRIPTION OF TARGETS 

Our first targeted S N, SN 2002bo, w as discov- 
ered on 2002 March 9 (jCacella et all l2002l ) and was 
spectroscop ically confirmed as a Type la SN shortly 
thereafter ('Kawa kita et al.l l2002t iBenetti et all 120021 : 

iMatheson et al 2002f ). The SN appeared approximately 

18" southeast of the nucleus of NGC 3190 in the mid- 
dle of an apparent dust lane. Early optical spectra of 
the SN indicated that it was quite young (10-14 days 
before maximum) and possibly reddened. Based on the 
SN's proximity to Earth and very early discovery, we 
triggered a pre-approved Chandra ToO observation on 
March 12. The nominally 20 ks long observation was 
carried out two days later (the midpoint of the observa- 
tion was at JD 2452347.9). We now know that the time of 
maximum light (in the B band) occurred o n 2002 March 
23.0 UT or, equivaleiitly JD 2452356.5, (jBenetti et all 
12004 iKrisciunas et al.l |2004[ ) and that the inferred ex- 
plosion epoch was JD 2452338.6 ± 0.5. Our Chandra 
ToO observation was therefore made only 9.3 days af- 
ter explosion. We adopt a value of 22 M pc for the dis- 
tance to NGC 3190 (for mo re details, see IBenetti et al.l 
12004 IKrisciunas et al.|[200l . SN 2002bo showed signif- 
icant reddening E{B -V) = 0.43 ± 0.10 (jBenetti et al.l 
|2004|) . most likely due to its host galaxy. We converted 
this to an equivale nt X-ray absorbin g column value of 
3 X 10^1 cm-2 using iCox et all (|2000l p. 527) in order to 
model the SN's X-ray flux. 

The second SN la we observed with Chandra^ 
SN 2002ic, is the first SN la in which a circumstcl- 
lar in teraction was detected (however, see Benctti et al. 
l2006l f or an alternate view) . Di scovered on 2002 Novem- 
ber 13 (jWood-Vasev et al]|2002| ). it was spectroscopically 
identified on 2002 December 7 as a t ype la SN at a 
redshift of 0666 fHa muv et aD 120021 ). In 2003 June 
iHamuv et al.l (2003a) reported the detection of broad 
(FWHM ^ 1800 km s~^) Ha emission from the super- 
nova, revealing a SN/CSM interaction in SN 2002ic in- 
volving a den se CSM with as much as several solar masses 
of material (jHamuv et al.l [2003b[ ). Evidence has been 
presented for a clumpy, aspherical st ructure for t he CSM 
llWang et al.ll200llDeng et al.ll200l FChugai & Ch evaheii 
I2007D as well as for gaps or troughs in the radial pro - 
file (|Wood-Vasev et al.ll200llChugai fc Yungelsonll200l . 
The unprecedented strength of the CSM interaction in 
SN 2002ic combined with the persistence of the Ha emis- 
sion (which lasted for many months), led us to request a 
20-ks observation with Chandra under the Director's dis- 
cretionary time program. The request was approved and 
our observation was carried out on JD 2452863.5 (the ob- 
servat ion began on 2003 August 11). IWood-Vasev et all 
(|2004l ) estimate a date for maximum light in the I3-band 
of JD 2452606. Assuming a typical ti me between explo- 
sion and maximum light of ^^18 days (jRiess "eFall [19991) 



we estimate an explosion date of JD 2452588. The Chan- 
dra observation was made 275 days after explosion or 
approximately 260 days in the SN frame. 

Another peculiar SN la (SN 2005gj), showing a circum- 
stellar interaction, was discovere d on 2005 September 
26 by the SD SS-II collaboration (jBarentine et al.ll2005l : 
lFriemanll2005f) . Originally classified as a probable type- 
Iln supernova because of the presence of resolved Hydro- 
gen line emission, further spectroscopy demonstrated a 
remarkable similarity with SN 2002ic, resulting in a re- 
classification of the event as a SN la but with clea r signs 
of a circumstellar interaction (|Prieto et al.ll200"5h . The 
SN lies in an anonym ous galaxy at a redshift of 0.062 
(jBarentine et al.ll2005f ). No X-ray emission was detected 
in a brief observation on 2005 November 24 using the X- 
ray t elescope onboard the Swift satellite (Imml er et al.l 
l2005f ). A 50-ks long Chandra observation was carried out 
under the Director's discretionary time (DDT) program 
on JD 2453 716.5 (the observatio n began on 2005 De- 
cember 11). lAldering et al.l (j2006t ) recently presented an 
analysis of the optical light curves and spectra in which 
they conservatively estimate that SN 2005gj exploded 
someti me in the interval September 18.6-24.6. More re- 
cently, iPrieto et al.l (j2007( ) estimate the time of explosion 
to be September 24.4 ± 2.0, based simply on whether or 
not they detected the SN, which, conservatively, can only 
provide an upper limit to the time of explosion. Since a 
difference of a few days is not signif icant for our pu r poses , 
we take the midpoint value of the lAldering et al.l (j2006l ) 
study (JD 2453635.1) as our estimate for the explosion 
date. This is some 77 days (in the SN frame) before the 
Chandra observation. 

The final SN la we report on here is SN 2005ke, which 
was discovered on 2005 November 13 (jPuckett et al.l 
120051) and confirmed as a Type la SN a few days later 
(Patat et al. 2005). These later authors estimated the 
spectral age by comparison to SN 1999by to be a few 
days before maximum light, which suggests an approxi- 
mate explosion date of Nov 2 (JD 2453676.50). Although 
the host galaxy, NGC 1371, is relatively nearby (we use 
a distance of 17 Mpc; iTullvl 1988) the SN was underlu- 
minous and not discovered particularly early so it was 
not initially observed by Chandra. However, Swift de- 
voted considerable observing time to the SN (~250 ks) 
over the course of the 2005/2006 winter and reported a 
tentative X-ray detection at about the 3 a significance 
level with a fltuc of ^4 x 10~ ^^ ergs cm~^ s~^ (0.3-2 
keV band) (jimmler et al.ll2006f l. The reported midpoint 
of the Swift observation is some 42 days after explosion. 
Subsequently a Chandra DDT observation was approved 
and carried out on JD 2453786.1 (beginning on 2006 
February 19) about 110 da ys after ex plosion. No radio 
emission has been detected (jSoderberg 2006). 

Table [1] summaries some key observational information 
about the several SNe discussed in this paper. 

3. OBSERVATIONS AND TECHNIQUES 

Each target was observed near the nominal telescope 
aimpoint on the back-side illuminated chip (S3) of the 
Chandra Advanced CCD Imaging Spectrometer (here- 
after ACIS-S). Although the background rate is higher 
here than on the front-side illuminated chips of the imag- 
ing array, the greatly enhanced low energy response of 
chip S3 was the main deciding factor in favor of using it. 



X-ray Flux Limits to Four Type la SNe 



3 



TABLE 1 

Observed Supernovae 





SN 2002bo 


SN 2002ic 


SN 2005gj 


SN 2005ke 


R.A. (J2000) 


10:18:06.5 


01:30:02.6 


03:01:12.0 


03:35:04.4 


Decl. (J2000) 


21:49:41 


21:53:07 


-00:33:13.9 


-24:56:38.8 


Explosion date (JD) 


2452338.6 


2452588.0 


2453635.1 


2453676.5 


Chandra Obs. date (JD) 


2452347.9 


2452863.5 


2453716.5 


2453786.1 


Chandra exposure (ks) 


19.6 


17.3 


49.4 


14.8 


A''hi (atoms cm~^) 


2.1 X lO^" 


6.2 X lO^^^ 


7.1 X lO^o 


1.4 X lO^o 


D (Mpc) / 2 


22 


0.0666 


0.062 


17 



Observations were carried out in timed exposure mode 
with events recorded in FAINT (SN 2002bo: Observa- 
tion Identifier [ObsID] 2760; SN 2005kc: ObsID 7277) or 
VFAINT (SN 2002ic: ObsID 4449; SN 2005gj: ObsID 
7241) data mode. The same basic data reduction pro- 
cesses were applied to each observation, using standard 
Chandra software and calibration (CIAO version 3.2.2, 
CALDB 3.1.0). This included applying the latest gain 
map, identifying and removing hot pixel and afterglow 
events, using VFAINT mode information (when avail- 
able) to remove some non- X-ray (i.e., charged particle) 
background events, and filtering on grade (retaining the 
usual values 02346) and status. Light curves of the en- 
tire S3 chip were made to identify times of high back- 
ground. The SN 2002bo, SN 2005gj, and SN 2005ke 
data were mostly free of flares, while the SN 2002ic data 
showed a modest amount of background flaring. Exclud- 
ing time intervals more than 3 a from the mean resulted 
in final useful exposure times of 19620 s (SN 2002bo), 
17280 s (SN 2002ic), 49410 s (SN 2005gj), and 14810 
s (SN 2005ke). The astrometry of each pointing was 
checked as well. We found no unambiguous matches 
between X-ray and optical point sources in the field of 
SN 2002bo, so it was not possible to obtain an inde- 
pendent check of the absolute astrometry of the X-ray 
data. In lieu of this the on-line "Aspect Calculator" 
tool was run, yielding an aspect offset of 0.5" which was 
applied to the data. For SN 2002ic, it was possible to 
match four X-ray sources with optical counterparts in the 
USNO-A2.0 catalog from which an aspect offset of 0.6" 
was determined for the Chandra data. There were nine 
optical counterparts to X-ray sources in the SN 2005gj 
field which resulted in an aspect correction of 1.0" to 
the Chandra data. One X-ray source showed a clear op- 
tical counterpart for the SN 2005ke observation. The 
positional offset was less than 0.1" so no adjustment was 
applied to the Chandra astrometry. 

Grayscale images of the Chandra data are presented in 
Figures [T]-[3] in two spectral bands: 0.5-2 keV (left) and 
2-6 keV (right). The small insert in the lower right corner 
of each panel shows the raw Chandra data in the imme- 
diate vicinity of each SN. There are no detected X-ray 
photons in a 1" (radius) circle centered on the positions 
of either SN 2002bo or SN 2005ke. Only two photons 
(with energies of approximate 1.1 keV and 3.6 keV) were 
detected at the position of SN 2002ic. Two photons (with 
energies of 4.0 keV and 7.7 keV) were also detected at 
the position of SN 2002gj These are clearly not significant 
detections and so in the following we demonstrate how 
we obtained upper limits to the X-ray flux of each SN. 
Since the range of allowed spectral models is large, we 
developed a technique that allows us to determine upper 



limits for just about any spectral form imaginable. 

One essential key to our method is to obtain an accu- 
rate spectral model for the background. To this end we 
extracted spectra from the large circular regions shown 
on Figures [THll excluding all serendipitous point sources 
and for each of SN 2002bo and SN 2002ke, an ellipti- 
cal region centered on the host galaxy that encompasses 
all the galaxy's apparent X-ray emission. The detector 
spectral response function (calculated using mkacisrmf) 
as well as the overall effective area function were deter- 
mined at a large number of detector positions within the 
extraction region and then weighted by the spatial dis- 
tribution of detected events to produce final weighted 
response functions. The black data points near the bot- 
tom of Figures [5] and [H] show the background spectra 
used for SN 2002bo and SN 2005ke (the other cases are 
similar to these two). Note that the spectra were re- 
binned to a signal-to-noise ratio of ~3 for display only; 
spectral fits were carried out on the original unbinned 
spectral data. The curve through these points indicates 
our best fit model, which consists of two power law com- 
ponents to describe the non- X-ray (i.e., particle) back- 
ground, three narrow lines to account for instrumental 
fluorescence lines from Si and Au, and two astrophys- 
ical models: a low temperature thermal model for the 
local Galactic emission and a hard power law for the un- 
resolved cosmic X-ray background. These later models 
included interstellar abs orptio n using the Galactic val- 
ues (Dickev & Lockman |l993) : TVh = 2.1 x cmr'^ 
(SN 2002bo), TVh = 6.2 x lO^^cm'^ (SN 2002ic), iVn = 
7.1 X 10^0 cm-2 (SN 2005gj), and A^h = 1-4 x W° cm-^ 
(SN 2005ke). One of the non- X-ray background power 
law models is required to model the spectrum below '^5 
keV; the other one fits the sharp upturn in emission be- 
yond 7 keV. Overall the non- X-ray background domi- 
nates the spectrum above about 2 keV. The fits were 
carried out using a maximum likelihood figure-of-merit 
function for Poisson-distributed data (the so-called "c- 
stat" option in xspcc). The background models for all 
four targets are similar as are the quality of the fits, 
which are acceptable from a statistical point-of-view (de- 
termined by examining the values on grouped data). 

Detector spectral response and effective area functions 
were calculated as above for the SN regions. We assume 
that the background intensity is uniform on arcminute 
spatial scales near the detector aimpoint and that its 
spectral form does not vary with position either. Thus 
we can scale the background models to the SN regions 
using the ratio of detector pixels in the source and back- 
ground spectral extraction regions. When this is done, 
we find that the background model predicts the detec- 
tion of 0.19 events (SN 2002bo), 0.20 events (SN 2002ic), 
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SN2002bO Chandra ACIS-S (0.5-2 keV) SN2002bO CAa/irfra ACIS-S (2-6 keV) 




R.A. (J2000) R.A. (J2000) 

Fig. 1. — Chandra ACIS-S X-ray images of the vicinity of SN 2002bo in the 0.5-2 keV (left) and 2-6 keV (right) bands displayed with 
1.968" pixel binning (block 4). The various circles and ellipses denote spectral extraction regions (see text). The small insert panel at the 
bottom right shows the raw data in the immediate vicinity of the SN in unblocked pixels (i.e., 0.492" X 0.492"). This image is centered 
on the SN position and the axis labels are in units of arcseconds as indicated. 



SN2002iC CAgnrfra ACIS-S (0.5-2 keV) 
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Fig. 2. — Chandra ACIS-S X-ray images of the vicinity of SN 
Otherwise similar to figure [T] 

0.50 events (SN 2005gj), and 0.14 events (SN 2005ke) 
over the entire Chandra band. The number of events 
actually detected in the observations of SN 2002bo and 
SN 2005ke (0 for each) is consistent with their respec- 
tive background rates. The Poisson probability to detect 

2 (or more) X-ray photons when 0.2 are expected, as is 
the case for SN 2002ic, is only 1.7%. The Poisson prob- 
ability to detect 2 (or more) X-ray photons purely from 
background in the case of SN 2005gj is 9%. 

We determine flux upper limits for a specific assumed 
astrophysical spectral model (including the measured 
redshifts for SN 2002ic and SN 2005gj), by adding it to 
the background model and then increasing its normaliza- 
tion until the c-stat value increases by A = 9, which is the 

3 a confidence level criterion. Results for two X-ray flux 
bands are shown in Tables [2 to [5] and plots of the varia- 
tion of the c-stat value vs. the total 0.5-10 keV band X- 
ray flux are shown in Fig. [T] All quoted fluxes have been 
corrected to remove Galactic absorption for SN 2002ic, 
SN 2005gj, and SN 2005ke, and the absorption from the 



SN2002iC Chandra ACIS-S (2-6 keV) 
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2005ic. The two circles denote spectral extraction regions (see text). 

host galaxy for SN 2002bo. The first set of four mod- 
els in these tables (dashed curves in the figure) assume 
thermal plasma emission with no intrinsic absorption, 
normal solar composition and temperatures from 2 to 
20 keV. The second set of four models (solid curves in 
the figure) assume a hot {kT = 80 keV) bremsstrahlung 
model with varying amounts of intrinsic absorption: from 
iVn = 1 X 10^2 cm~2 to TVh = 8 x 10^^ cm'^. Uncertainty 
in the level of the background has a very small effect on 
the flux upper limits. For example, assuming a 20% vari- 
ation in the level of the model background for SN 2002ic 
introduces a ~3% change in the derived flux limit for the 
hard spectral model. 

Fig. [8] graphically illustrates how the upper limits de- 
pend on the assumed spectral shape. Each curve here 
corresponds to the upper limit for that spectrum allowed 
by the Chandra data. Clearly more hard band emission 
is allowed if the soft band is depressed, and vice versa. 

Some comments on Fig. [7] are warranted. The top 
left panel (SN 2002bo) shows a steady rise in the c-stat 
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Chandra ACIS-S X-ray images of the vicinity of SN 2005gj. Otherwise similar to figure [T] 



SN2005ke 



Chandra ACIS-S (0.5-2 keV) 



SN2005ke 




Chandra ACIS-S (2-6 keV) 
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Fig. 4. — Chandra ACIS-S X-ray images of the vicinity of SN 2005ke. The various circles and ellipses denote spectral extraction regions 
(see text). Otherwise similar to figure [T] 
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Fig. 5. — ACIS-S spectra of the background (black, bottom) 
and several regions in NGC 3190: outer galaxy (blue), mid galaxy 
(green), inner galaxy (red), and nucleus (black, top). 
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Fig. 6. — ACIS-S spectra of the background (black, bottom) and 
several regions in NGC 1371: outer galaxy (blue), inner galaxy 
(red), and nucleus (black, top). 
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10-16 10-16 10-14 10-16 10-1' 

X-ray flux (ergs cm^^ s"') X-ray flux (ergs cm^^ s"') 



Fig. 7.— X-ray flux limits in the 0.5-10 keV band for SN 2002bo (top left), SN 2002ic (top right), SN 2005gj (bottom left) and 
SN 2005ke (bottom right). The likelihood fit statistic ("c-stat") is plotted as a function of the X-ray flux limit. Solid curves come from 
assuming thermal plasma models with different temperatures: kT = 2 (red), 5 (blue), 10 (green), and 20 (black) keV. Dashed curves come 
from assuming a hot (kT = 80 keV) bremsstrahlung emission model with different amounts of intrinsic absorption: TVjj = 1 X 10^^ (red), 
A^H = 2 X 10^2 (blue), A^H = 4 X 10^^ (green), and A^h = 8 X 10^^ (black) atoms cm ■2. The horizontal dotted line near the top of each 
panel denotes the S-c criterion on the likelihood fit statistic. The intersection of this line with the appropriate curve provides the 3-(t X-ray 
flux upper limit (to be read off the x-axis). 



TABLE 2 

X-RAY Flux Limits (3 a) for SN 2002bo 



kT 


Intrinsic Afjj 


Flux 


Flux 


(keV) 


(atoms cm^2) 


(0.5-2 keV) 


(2-10 keV) 


2 




1.2 X 10-15 


6.8 X lO-i'' 


5 




9.9 X IQ-^^ 


1.7 X 10-15 


10 




9.0 X 10-1"^ 


2.2 X 10-15 


20 




8.5 X IQ-^^ 


2.5 X 10-15 


80 


1 X 10^2 


3.2 X lO-i** 


4.1 X 10-15 


80 


2 X 10^2 


2.0 X lO-i** 


5.4 X 10-15 


80 


4 X 1022 


8.3 X lO-i'^ 


7.1 X 10-15 


80 


8 X 1022 


1.4 X lO-i'^ 


9.2 X 10-15 



value as the flux (or normalization) of the assumed spec- 
tral model is increased. The top right panel (SN 2002ic) 
shows a slightly different behavior: the c-stat value ini- 
tially decreases as the flux of the model increases until it 
reaches a minimum after which it begins to grow rapidly. 
What is happening is that the added spectral model is 
accounting for the 2 detected photons in the SN 2002ic 
spectrum. The reduction in the c-stat value is A ~ 8 
which is significant at between 2 and 3 a. The curves for 
SN 2005gj resemble those for SN 2002ic for essentially the 



TABLE 3 

X-RAY Flux Limits (3 a) for SN 2002ic 



kT 


Intrinsic Nn 


Flux 


Flux 


(keV) 


(atoms cm— 2) 


(0.5-2 keV) 


(2-10 keV) 


2 




3.0 X 10-15 


1.5 X 10-15 


5 




2.7 X 10-15 


4.2 X 10-15 


10 




2.5 X 10-15 


5.8 X 10-15 


20 




2.4 X 10-15 


6.9 X 10-15 


80 


1 X 1022 


1.1 X 10-15 


1.5 X lO-i"! 


80 


2 X 1022 


6.5 X 10-1^ 


1.7 X lO-i"! 


80 


4 X 1022 


2.4 X 10-16 


2.0 X 10-14 


80 


8 X 1022 


3.8 X 10-" 


2.5 X lO-i"! 


Chugai 


Model - Day 260 


9.0 X 10-1^ 


2.2 X 10-1" 



same reason, while the curves for SN 2005ke (for which 
no photons were detected) resemble SN 2002bo. 

Finally we compare our X-ray flux limi ts for SN 2005gj 
with those quoted bv lPrieto etall ()2007f ): 5 x IQ-^^ ergs 
cm-^ s"^ (at 68% confidence level) and 9 x 10"^^ ergs 
cm-^ s-^ (at 95.5% confidence level) for the 0.5-8 keV 
band. The spectral form they assumed (an unabsorbed 
power law with photon index T = —2) is roughly compa- 
rable to our thermal plasma model with kT — 10 keV. 
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TABLE 4 

X-RAY Flux Limits (3 a) for, SN 2005gj 



kT 


Intrinsic Afjj 


Flux 


Flux 


(kcV) 


(atoms cm^^) 


(0.5-2 kcV) 


(2-10 kcV) 


2 
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TABLE 5 

X-RAY Flux Limits (3 a) for SN 2005ke 



kT Intrinsic A^h 
(keV) (atoms cm-2) 



Flux 
(0.5-2 keV) 



Flux 
(2-10 keV) 



2 




9.4 X 10- 


16 


5.4 X 10- 


16 


5 




8.5 X 10- 


16 


1.4 X 10- 


15 
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16 
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17 
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Fig. 8. — Maximum allowed X-ray emission (at 3 a) from 
SN 2002ic for several different assumed spectral forms. Curves with 
emission lines are for thermal models with normal (solar) composi- 
tion, no intrinsic absorption, and a range of different temperatures: 
kT = 2 keV (cyan), kT = 5 keV (red), kT = 10 keV (green), 
kT = 20 keV (blue). Smooth curves are for a hard bremsstrahlung 
model {kT = 80 keV) and varying amounts of intrinsic absorp- 
tion: ATjj = 1 X 1022 cm-2 (cyan), A^h = 2 X 1022 cm-2 (red). 
Ah = 4 X 1022 cm-2 (green), and A^h = 8 X 1022 cm-2 (tilue). 

We utilize Fig. [7] to convert to the 68% (c-stat A = 1) 
and 95.5% (c-stat A = 4) confidence levels and obtain 
comparable flux limits (0.5-10 kcV band) of 4 x 10 

'■^ and 9 x 10^^^ ergs cm-^ s~^, respectively. 
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3.1. Re-evaluation of tentative Swift X-ray detection of 
SN 2005ke 



The importance of the tentative detection of SN 2005ke 
by the Swift X-ray Telescope (XRT) ins trument can- 
not b e overstated, yet the original report (|Immler et al.l 
l2006f ) provided only a brief summary of the data reduc- 
tion and analysis. We undertook an independent analy- 
sis of these data in order to address some issues that we 
considered to be incompletely resolved in the published 
study. 

We downloaded the archival Swift XRT data on 
SN 2005ke observed during the winter of 2005/2006 
from the HE AS ARC. There were 43 separate obser- 
vations from 14 November 2005 until 5 March 2006 
encompassing ObsIDs from 0030341001 to 0030341008 
and from 0030342001 to 0030342036. Recently (March 
2007) Swift reobserved SN 2005ke (ObsIDs 0035898002 
to 0035898004) and we examined these data as well. The 
report on this recent observation can be summarized suc- 
cinctly: no photons were detected within 14" of the SN 
position. In our analysis of the rest of the data we began 
with the standard cleaned events files utilizing only data 
in photon counting mode and all good grades (0-12). 

The position of the nuclear X-ray source was located 
in each of the separate observations and then images in a 
broad band (0.5-6 keV, pi bins 50:600) were generated us- 
ing this location for the image center. The four shortest 
observations (exposures < 100 s) gave no detected pho- 
tons from the nucleus and were excluded from further 
consideration. The total dead-time corrected exposure 
in the merged image was 251.0 ks. 

We registered the Swift data to the Chandra data using 
the observed positions of 11 fairly bright X-ray sources 
that contained at least 30 counts in the merged Swift 
image. A simple linear 4" shift in position was sufficient 
to bring the Swift and Chandra positions into agreement 
with a scatter in the relative after-shift positions of <1". 
The accuracy of the Chandra positional astrometry was 
verified, as mentioned above, using the cataloged posi- 
tion of an optical star detected in the X-rays. Figure [5] 
shows the Swift (left panel) and Chandra (middle panel) 
broad band images highlighting the region covering the 
supernova position and the nucleus of NGC 1371. 

The broadband Swift and Chandra count rates were 
determined for the aforementioned 11 sources in addi- 
tion to 9 other fainter ones. These are essentially all the 
compact X-ray sources within ~5' of the Chandra opti- 
cal axis. For Swift the rates were extracted from within 
9" radius circular regions and then multiplied by two 
to account for the point-spread function (PSF) of the 
telescope (the half-power diameter at 1.5 keV is 18"). 
The Chandra point source extraction regions were large 
enough to capture >90% of the source flux. Figure [10] 
plots the Swift vs. Chandra count rates. The rates are 
well correlated and within the errors most of the sources 
(14/20) are consistent with the best-fit average ratio of 
rates (Swift to Chandra) of ~0.19 (shown as the long 
dashed line in the figure). This is close to the ratio of 
effective areas at 1.5 keV for the two observatories (135 
cm^ / 600 cm^ k, 0.225), which is shown as the short 
dashed line in the figure. The six sources that deviate 
from the average ratio of rates do so by generally no 
more than about a factor of two and are likely just time 
variable sources. 

Having demonstrated the astrometric and photomet- 
ric accuracy of the Swift data, we are now in the po- 



8 



Hughes et al. 



SN2005ke 



Swift XRT (0.5-6 keV) 



SN2005ke 



Chandra ACIS-S (0.5-2 keV) 



NGC 1371 



H alpha 



R.A. (J2000) 



R.A. (J2000) 



R.A. (J2000) 



Fig. 9. — The left and middle panels show X-ray images of the vicinity of SN 2005ke from the Swift XRT (left) and Chandra ACIS-S 
(middle). The solid-linetype circle near the middle of these panels denote the extraction regions for SN 2005ke (9" radius for Swift, 1" radius 
for Chandra). In the case of Swift the small dotted-linetype circles denote the locations from which background was obtained. Note that the 
4 background and single source regions are equidistant from the center of the bright nuclear emission (as note d by the large dashed-l i netyp e 
circle which has a 1' radius). The right panel shows the Hq emission from NCG 1371 originally presented in lHameed fc Devereuxl 1)19991 ). 
All three images cover the same field of view. 

ceeds the background by more than 2 a. This is not a 
/ significant detection, so we determine a count rate up- 

^ per Hmit of 1.6 x 10~^ s~^ (at 3 sigma confidence) from 

: the Swift data. This corresponds roughly to a Chandra 

count rate of 8 x 10~^ s~^ (using a ratio of rates of 0.2), 
which is several factors higher than the Chandra count 
'77 rate upper limit we derived for SN 2002bo. 



0.001 



rf, 0.0001 



0.5-6 keV band 



0.0001 



0.001 

Chandra rate (s"^) 



0.01 



Fig. 10. — Chandra ACIS-S vs. Swift XRT count rates and 1 a 
error bars for unresolved X-ray sources in the field of view of the 
SN 2005ke observations. The long dashed line shows the mean ratio 
of rates from the measured data points. This is in good agreement 
with the ratio of effective areas at 1.5 keV, which is shown as the 
short dashed line. The hatched region across the middle of the 
plot shows the clai med detection count rate and 1 a error from 
llmmler et al.l II2006I) , while the thick dark line shows our 3 a upper 
limit to the count rate. 

sition to determine the X-ray count rate and detection 
significance of SN 2005ke. First w e note that th e count 
rate of 2.2 x 10""* s^^ quoted by llmmler et all ()2006D 
would lead us to expect something like 25 detected pho- 
tons, after background subtraction, within our fiducial 
9" radius circular region. In fact, we detect only 15 total 
events, which includes both source and background emis- 
sion. We estimate the background using the 4 regions de- 
noted with dotted circles in Fig.[9l which are at least 18" 
from the SN position, are all equidistant ('^1') from the 
center of the nuclear emission (as indicated by the large 
dashed circle) and have the same radius as the source 
region. This yields an estimate of 7.8 ± 1.4 background 
events for the SN 2005ke source region. Our background 
subtracted count rate (including the multiplicative PSF- 
correction factor of 2) becomes (5.7 ± 3.3) x 10^^ s^^. 
Even if we vary the source extraction region (by doubling 
it, for example), the number of source counts never ex- 



4. LIMITS ON THE AMBIENT CSM 



We use the model of IChugai et al.l ()2004l . hereafter 
CCL) to determine constraints on the density of the am- 
bient medium surrounding SN 2002bo, SN 2002ic, and 
SN 2005gj from our X-ray flux upper limits. We do not 
consider SN 2005ke, since the time of X-ray observation 
(>40 days after explosion) and the limit on X-ray lumi- 
nosity are both less constraining of the SN models than 
is the case for SN 2002bo. Here we provide a brief sum- 
mary of the calculations; the reader is referred to CCL 
for more details. 

The model considers interaction of the SN ejecta with 
spherically symmetric, smooth circumstellar gas and 
treats the dynamics in a thin shell approximation. The 
evolution is followed by solving the equations of motion 
numerically assuming an exponential density profile for 
the SN envelope and a power law CSM density profile 
(which in the case of SN 2002ic includes different indices 
for different radial ranges). Specifically, the density of 
SN la ejecta is described by the law p = po exp(— w/wq), 
where po and vq are defined by the mass M and kinetic 
energy E. We adopt a typical ejecta mass M = 1.4 Mq, 
an d kinetic energy E = 1-4 x 10^^ ergs (the model PDD3 
of iHoflich &: Khokhlovl 11996^ . The outermost unburnt 
layers of the ejecta {v > 15000 km s~^) are represented 
as pure oxygen. Interior to this the ejecta are Si- and Fe- 
rich, while the CS wind is assumed to be hydrogen-rich. 
The numerical model provides the time evolution of the 
shock radii, velocities, and temperatures. 

The internal structures of the forward and reverse 
shocks are not determined in the model, so the X- 
ray luminosity at age t is calculated from the instan- 
taneous kinetic luminosity of the shock, = T^Lkin 
(jChevalier fc FranssonI 11994 IChugail I1992D . The radia- 
tion efficiency factor rj is given hy r] = t/{t + tc), where tc 
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is the cooling time of the post-shock gas. Electron and 
ion temperatures are fully equilibrated for SN 2002ic and 
SN 2005gj due to Coulomb collisions (however, see below 
for SN 2002bo) and the shape of the emission is given 
simply by a thermal bremsstrahlung spectrum at the 
appropriate temperature; no line emission is included. 
This is a conservative approximation, since including line 
emission would boost the modeled X-ray emission and 
thereby decrease our limit on the ambient gas density. 

The reverse shock is essentially radiative in all cases 
we consider here, so the shocked metal-rich SN ejecta 
create a cool dense shell (CDS) at the contact disconti- 
nuity, which causes severe absorption of X-ray emission 
from the reverse shock. The absorption coefficient for the 
metal-rich ejecta is given by kx ~ 5000(i?/l keV)"®/'^ 
cm^ gm^^. X-rays are also absorbed by the CSM, al- 
though this material is assumed to be hydrogen-rich for 
which kx is some 50 times smaller. That portion of the 
total X-ray luminosity absorbed by SN unshockcd ejecta, 
CDS, and unshocked circumstellar (CS) gas is presumed 
to be fully re-emitted in the optical band. 

One of the most important quantities influencing the 
X-ray luminosity of the SN/CSM interaction is the den- 
sity of the wind. Typically this has been expressed in 
terms of a constant wind speed, fm, and mass loss rate, 
M, as pw = M / (Anvwr'^) ■ The quantity that we con- 
strain with the Chandra observations is the wind density 
parameter w = M When quoting mass loss rates we 
will give values in terms of ti^.io = Uu,/10kms~^. 

4.1. SN 2002ho 

It is doubtful that the shocked electrons and ions in the 
forward shock reach temperature equilibration in the low 
density CS wind around SN 2002bo. In general the av- 
erage post-shock electron temperature, Te, will be in the 
range /ie//J < T^/Ti < 1, where /ie — nie/rnp — 1/1840 
and pL is the mean molecular weight. We calculate in- 
termediate values of the electron temperature using the 
relation Te/T, = 1 - [1 - {pLe/fi)][l - exp(-6l)]/6l, where 
9 = t/toq and toq, the Coulomb equilibration time, is 
evaluated conservatively assuming Tg = T^. This expres- 
sion is the average T^/Ti for an adiabatic plane shock 
in a homogeneous medium at age t for a constant post- 
shock equilibration timescale. Because we use the max- 
imal value of tcq, the value of recovered in this way 
will be underestimated, and accordingly, our wind den- 
sity constraints will be somewhat overestimated. 

We ran a set of models with wind density parame- 
ter values in the range 6 X lO^-* g cm-i to 3 X lO^^ 
g cm~^ without assuming complete temperature equili- 
bration for a SN age o f 9.3 d and distance of 22 Mpc 
(|Krisciunas et al.l l2004j) . Each trial spectrum was in- 
put to the spectral fitting program and compared to the 
Chandra data in exactly the same manner as described 
above (§|31)- A value oiw = 1.2 x 10^^ g cm~^ was found 
to be the maximum value allowed by the Chandra data 
at the 99% confidence level. Fig. [11] plots this model 
along with a rejected case with higher wind density pa- 
rameter and, for reference, we plot the broadband flux 
upper limits for the absorbed bremsstrahlung model with 
Nh = 2x 10^2 atoms cm'^. 

At about 10 keV the model flux is dominated by the 
reverse shock emission, while the forward shock with its 
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Fig. 11. — Model X-ray emission from SN 2002bo on day 
9.3 after explosion. The upper panel shows the case with non- 
equilibrated temperatures (Te < Ti) for wind density parameters 
of to = 1.2 X 10^^ g cm~^ (solid curve) and «) = 3 X 10^^ g cm~^ 
(dashed curve), corresponding, respectively, to models "neql2" and 
"neqSO" in tal3le[6] The lower panel shows the model spectra with 
equal electron and ion temperatures for ui = 6 X 10^'' g cm~^ 
(solid curve, model "eqO.6") and w = 1.5 X 10^'* g cm~^ (dashed 
curve, model "eql.5"). In each panel the solid curve corresponds 
to the 99% upper limit on X-ray emission; the dashed curve indi- 
cates how the spectra vary with increasing wind density parameter. 
For comparison we plot representative Chandra upper limits in two 
broad energy bands using the absorbed bremsstrahlung model with 
Nh = 2 X 10^^ atoms cm"^. 

cooler electrons produces the weaker peak at about 1 
keV. The low electron temperature of the forward shock 
is the result of inefficient equilibration in the low density 
wind. The electron and ion temperatures in the reverse 
shock, on the contrary, are close to equilibrated. More- 
over, the reverse shock is radiative during the initial pe- 
riod up to '^28 days for wind parameters in the range 
being considered. The cool dense shell of the reverse 
shock is responsible for the strong absorption of X-rays 
from the reverse shock. 

Given that the amount of coUisionless electron heating 
at the forward shock is highly uncertain we also ran a 
set of models for the equipartition case, Tg = Ti. Here 
the maximal wind density parameter is w = 6 x 10^^ g 
cm~^ (Fis. [TT|l. Remarkably, in the latter case the peak 
at ~1 keV is dominated by Comptonization of the SN 
radiation by the hot electrons of the forward shock. The 
Comptonizcd spectral component is described by a power 
law oc E^°' with a depende nce on tt and fcTe/rrieC^ taken 
according to Pozdnvako v. Sobol. fc SunvaevI ()1976f ) . 

For reference in Table[S]we present numerical values for 
the parameters corresponding to the interaction models 
for SN 2002bo on day 9.3 that are plotted in Fig. [ill 
The top two models correspond to the non-equilibrated 
shocks {Te < Ti), while the other two models correspond 
to the equilibrated ones (Te = Ti). The columns list: a 
model designation, the wind density parameter (w), thin 
shell radius (i?s), velocity of the reverse (urs) and for- 
ward (t^Fs) shocks, electron temperature of the reverse 
(?e,Rs) and forward (Te.Fs) shocks, hydrogen column 
density of the CS gas (iVn) for r > Rs assuming H abun- 
dance X = 0.7. The last column gives the mass of the 
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cool dense shell formed by the shocked SN ejecta. This 
metal-rich matter strongly suppresses X-rays from the 
reverse shock and partially absorbs X-rays from the rear 
side of the forward shock. 

We also applied Lundqvist's models for the X-ray emis- 
sion emerging from the presumed wind surrounding a SN 
la (Lundqvist et al. 2007, in preparation). Note that 
there are some differences in assumptions and details 
in this model compared to those just presented. How- 
ever our intent here is not to resolve these differences, 
but rather to present an independent estimate of the ex- 
pected X-ray flux from the interaction between SN la 
ejecta and CSM at e arly times. Lund qvist's models use 
a similarity solution (jChevalierj 1 1982ar ) for the evolution 
with power law indices of n = 7 in the ejecta and s = 2 
in the wind. The velocities of the forward and reverse 
shock at 9.3 days after explosion are 2.9 x lO'' km s~^ 
and 7.0 x 10^ km s~^, respectively. The full ejecta den- 
sity and temperature structure from the similarity solu- 
tion is used for calculating the X-ray emission, assuming 
that the composition of the outer ejecta is equal parts 
of carbon and oxygen. We verified that the density in 
the reverse shocked ejecta is sufficient to assure that the 
electron and ion temperatures are equilibrated there. In 
these models a radiative cold dense shell does not form 
and so it is emission from the reverse shock that domi- 
nates (and in fact these calculations do not include for- 
ward shock emission). Only thermal emission is included; 
although inverse Compton scatteri ng could be an impor- 
tant contribution (see, e.g.. Fig. 7 in'Fransson et al.]ll996L 
for the case of SN1993J), we ignore it here (which is again 
a conservative assumption as regards the limit on CSM 
density) . Intrinsic absorption comes about only from the 
overlying CSM. Models were generated and compared to 
the Chandra spectrum as above. The wind density limit 
we obtain is comparable to those above: w — 1.6 x 10^^ 
g cm~^. 

4.2. SN 2002ic 

SN 2002ic has been studied extensively. We take the 
position here that the intense late-time {t > 50 d) opti- 
cal luminosity of SN 2002ic is powered by the interaction 
of SN la ejecta with a dense circumstella r environment 
(iHamuv et al.l l2003bl: IWang et all l2004t IChugai et all 
120041 iNomoto et al.ll2005f ). The CCL interaction model 
that accounts for the optical luminosity and spectrum re- 
quires a dense stellar-wind-type density profile suggest- 
ing a high mass-loss ra te (M ~ 10~^i',i,in Mg yr~M from 
the progenitor system (jChugai fc Yungelsonll2004l ) . This 
model predicts an X-ray flux of order 10~^ photon cm~^ 
s^^ kev~^ in the 3-10 keV energy band for ages in the 
range 200 — 400 d. The predicted X-ray flux on day 
260 along with contributions of the reverse and forward 
shocks for the standard parameter set (CCL) is plotted 
in Fig. [12] for a distance 285 Mpc along with an esti- 
mate of the Chandra flux upper limit (using the most 
absorbed bremsstrahlung model with Nh = 8 x 10^^ 
atoms cm~^). We dub this model 'unmixed' for reasons 
that will become clear below. The modeled X-ray flux is 
clearly dominated by the emission of the forward shock; 
emission from the reverse shock is absorbed by the cool 
dense shell formed in the radiative ejecta. We compared 
this spectrum to the Chandra data and found that the 
X-ray flux exceeds the observational upper limit by a 
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Fig. 12. — The predicted emergent X-ray spectrum from 
SN 2002ic on day 260. For comparison we plot representative 
Chandra upper limits in two broad energy bands using the most 
absorbed bremsstrahlung model with N[j = 8 X 10^^ atoms cm~^. 
The upper panel shows the total flux (thick solid), flux from the 
forward shock (dotted), and flux from the reverse shock (thin solid) 
predicted by the model of CCL (see text). This model overpredicts 
the X-ray emission by a factor of ~4.5. The lower panel displays 
the total X-ray flux computed for the same model as above but with 
fragments of cool dense ejecta homogeneously mixed in the forward 
shock region for different values of the occultation optical depth: 
Toe = 1 (dashed line), Toc = 2 (thick solid line), Toc = 5 (thin 
solid line), and Toc = 10 (bottommost thin solid line). The dashed 
curve is rejected by the Chandra observation and the Toc ~ 2 case 
corresponds to the 99% allowed upper limit. 



factor of « 4.5. The disagreement is significant and indi- 
cates that the model is incorrectly predicting the X-ray 
emission from the forward shock. 

The model ignores possible deviation fro m spherical 
symm etry suggested by polarization data (|Wang et al.l 
l2004f ). Generally, if the CS gas is concentrated in a dense 
thin equatorial disk (thickness <^ radius), it is conceiv- 
able that the bulk of the kinetic energy of the interaction 
could be dissipated deep inside the SN envelope. As a 
result, the ratio of the emergent X-ray flux to the opti- 
cal flux could be small because of strong absorption in 
the dense disk. However this is unlikely to be the case. 
Indeed, large deviations from global spherical symme- 
try for the interaction are ruled out by the similarity of 
line widths and proflles for all known SN la with strong 
CS interaction, i.e., 2002ic, 1999E, 1997cy (CCL). Fur- 
thermore a modest level of asymmetry should not signifl- 
cantly affect the ratio of the emergent-to-absorbed X-ray 
luminosity compared to the spherical case. 

More serious consequences for the X-ray luminosity 
might result from the omission of the internal structure 
of the forward shock in the thin shell approximation. 
To estimate the error introduced by this approximation 
we computed a hydrodynamical model for the standard 
parameter set as before (CCL) using a one-dimension 
hydro dynamical Lagrangian code with artifical viscosity 
(e.g.. lGull|[l973l ). On day 260 the radius of the con- 
tact discontinuity of the hydrodynamical model coincides 
with the radius of the thin shell model to an accuracy of 
1%. The forward shock structure of the hydrodynami- 
cal model was used to calculate the X-ray emission from 



X-ray Flux Limits to Four Type la SNe 



11 



TABLE 6 
SN 2002BO Model Parameters 



Model W R, V^s ^FS Te,RS Te,FS A^H McDS 

(lO^gcm^^) (10^^ cm) (10® cm s'^) (10* cm s^^) (keV) (keV) (10^^ cm'^) (10"^ M©) 

neql2 12 2.22 2.68 24.9 18.3 1.93 18 5.6 

neqSO 30 2.02 2.72 22.4 22.6 4.76 50 11 

eqO.6 0.6 2.81 2.18 32.7 16.5 1290 0.71 0.33 

eql.5 1.5 2.64 2.41 30.4 20.1 1110 1.9 0.92 



the shocked gas assuming instant equihbration of ion and 
electron temperatures. The X-ray luminosity of the for- 
ward and reverse shocks is lower by factors of 1.5 and 
1.3, respectively, compared to the thin shell model, while 
the temperature of the forward shock is greater (80 keV 
compared to 50 keV in the thin shell model). The lower 
X-ray luminosity of both shocks in the hydro model re- 
sults in an optical luminosity lower by a factor of 1.3. To 
recover the optical bolometric luminosity the CS density 
in the hydro model needs to be slightly larger. We found 
that a R:!20% increase in the CS density is sufficient to 
recover the required optical luminosity. The new total 
mass of the CS envelope is Ril.9 Mq within a 3 x 10^^ 
cm radius compared to «1.6 Mq found in the thin shell 
approximation (CCL) . The increase of the CS density by 
20% slightly increases the X-absorption thus reducing the 
peak of the X-ray flux by less than 20%. We thus con- 
clude that differences between the hydrodynamical and 
thin shell models are rather minor and cannot account 
for the significant disparity between the Chandra upper 
limit and the predicted X-ray flux. 

Two missing factors that could potentially resolve 
the issue are the Ray leigh- Taylor (RT) instability of 
the CDS and clumping of the CS gas. The RT in- 
stability of the shocked decelerated ejecta is a generi c 
attribute of the SN/CS interaction (jChevahed ri982bf ). 
while the dumpiness of the CS matter around SN 2002ic 
has been invoked already to account for the Ha line 
profile (CCL). The RT instability leads to fragmen- 
tation and mixing of CD S fragments with the hot 
gas of the forward sho ck ('Chevali er fc BlondinI Il995l : 
iBlondin &: EllisonI [200l[ ) . Clmnpincss of the CS matter 
favors more efficient mixing of the CDS fragments within 
the forward shock principally through two mechanisms. 
First, the CS clouds penetrate deep inside the intercloud 
shock before they get shocked , fragmented a nd mixed 
(|Klein. McKee. fc Colenal[l99l lBlondinl[2?)0l . As a re- 
sult the bulk of the kinetic energy related to CS clouds 
dissipates closer to the contact discontinuity. Second, 
the CS clouds engulfed by the intercloud shock gener- 
ate vortex turbulence that favors more efficie nt penetra- 
tion of the RT spikes in the forward shock pun at al.l 
We thus propose a scenario for SN 2002ic in which 
the forward shock is no longer a regular layer of hot gas 
but instead is a macroscopic mixture of hot gas and cool 
metal-rich fragments of the CDS. The major outcome of 
this modification should be an additional component of 
absorption of the X-ray emission from the forward shock 
by the intermixed fragments of metal-rich CDS. 

To illustrate this effect we will calculate the absorp- 
tion of X-rays from the forward and reverse shocks as- 
suming the same thin shell model as in CCL, but with 
the following corrections. The CDS is assumed to be 



fragmented and homogeneously mixed in the forward 
shock layer. Remarkably, the analysis of the optical spec- 
trum of SN 2002ic suggests that the CDS is indeed frag- 
mented and well mixed. The ratio of the cumulative 
area of clumps to the area of the spherical shell, which 
is a measure of mixing, was estimated to be S/ Sq ~ 50 
on day 234 (CCL). The occultation optical depth of the 
fragmented CDS, i.e., the average number of fragments 
on the line of sight in the shell for randomly oriented 
plane fragments is then Tqc = 0.5(5'/5'o) ~ 25. Small 
values of the occultation optical depth, Tqc 1 — 2, 
would mimic the possible presence of "holes" in the mix- 
ing zone due to strong angular var iation in the column 
density (jChevalier fc BlondinI I1995D . Another relevant 
value is the average optical depth for X-rays r, which 
is determined by the average column density of CDS 
material and the X-ray absorption coefficient, for which 
we take the metal-rich case as given above. The effec- 
tive optical depth of the ensemble of fragments is then 
Toff — Toc[l — exp(r/Toc)]. In the limit Tcfi ^ 1 the effect 
of mixing for a homogeneous thin spherical layer with 
unabsorbed luminosity Lq is to reduce the emergent X- 
ray luminosity to a value of Lq/ (Atcs). Indeed, only half 
of the photons emitted in a layer with thickness equal 
to the mean free path (AR/Tes) can escape. The other 
factor of one-half comes from averaging over the angle 
between the photon direction and surface normal. 

The mixing we invoke may, however, be incomplete, 
which can be mimicked by adopting a lower occultation 
optical depth of the homogeneous layer, Tqc < 25. We 
also take into account additional absorption by the frag- 
mented shocked CS clouds with normal composition, as- 
suming that 50% of the swept-up CS gas resides in the 
cool shocked cloud material with Tqccs = 1- This compo- 
nent, presumably responsible for the Ha emission, pro- 
duces negligible X-ray absorption. 

The computed total flux of escaping X-rays is pre- 
sented in the lower panel of Fig. [T^] for Tqc = 1, 2, 5 and 
10, which vary from low to high degrees of mixing of CDS 
fragments into the forward shock region. The plot shows 
a pronounced effect on the X-ray absorption (see Fig. [T^ 
by the macroscopically mixed CDS material. The case 
Toe — 2 corresponds to the maximal ffux tolerated by the 
Chandra upper limit. This study indicates that mixing 
of the shocked SN ejecta with the hot gas of the forward 
shock can resolve the X-ray flux disparity between the 
original model that explains the optical spectrum (CCL) 
and the new Chandra observations of SN 2002ic. 

4.3. SN 2005gj 

According to its early spectrum and photometry, 
SN 2005gj is remarkably similar to SN 2002ic and very 
likely is a member of this new class of bright la SNe 
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Fig. 13. — The model X-ray spectrum of SN 2005gj. For com- 
parison we plot representative Chandra upper limits in two broad 
energy bands using the highly absorbed bremsstrahlung spectrum 
with A'^^f = 8 X 10^^ atoms cm~^. The upper panel shows the 
unmixed model used previously for SN 2002ic and recomputed at 
an epoch of i ~ 80 d as appropriate for SN 2005gj. This model 
overpredicts the X-ray emission by a factor of ~4. The lower panel 
displays the total X-ray flux computed for the same model as above 
but with fragments of cool dense ejecta homogeneously mixed in 
the forward shock region for different values of the occultation 
optical depth; Toc = 1 {dashed line), Tqc = 2 {thick solid line), 
Toe = 5 {thin solid line), and Toc = 10 {bottommost thin solid 
line). The dashed curve is rejected by the Chandra observation 
and the Toc ~ 2 case corresponds to the 99% allowed upper limit. 

embedded in dense CS envelope. This position has been 
str engthened by the recen t hght curve and spectral study 
of lAldering et al.l ()2006D . Their V- and I-band Ught 
curves for SN 2005gj m atch the shape of the early time 
evolutionary models of iChugai fc YungelsonI (|2004.) for 
SN 2002ic, although the models need to be scaled upward 
in brightness by 0.5 mag or so. The CSM interaction 
in SN 2005gj is evidentally stronger than in SN 2002ic. 
In the absence of a detailed hydrodynamical model for 
SN 2005gj, we utilize the same model and parameter val- 
ues as for SN 2002ic. First we computed the X-ray flux at 
the Chandra observation epoch {t ~ 80 d) neglecting the 
effects of ejecta mixed into the forward shock. The result- 
ing spectrum is presented in the upper panel of Fig. fTSl for 
a distance to SN 2005gj of 266 Mpc together with an esti- 
mate of the Chandra broad band flux upper limits (using 
the most highly absorbed bremsstrahlung spectrum with 
Nh = 8 X lO^'^ atoms cm~^). The model spectrum over- 
predicts the X-ray flux by a factor of 4. Although little is 
cur rently published about the environment of SN 2005gj , 
the lAldering et al.l ()2006( ) study suggests that the den- 
sity of the ambient medium is, if anything, higher than 
that around SN 2002ic, which exacerbates the discrep- 
ancy with the Chandra flux limit. Therefore we next 
considered the mixing of CDS fragments into the forward 
shock region in the case of SN 2005gj. The required oc- 
cultation optical depth to not exceed the Chandra limit 
is Toc < 2 (Fig.[T31). Consequently, if the CS environment 
around SN 2005gj is as dense as that around SN 2002ic, 
mixing of the fragmented CDS must be already signifi- 
cant at this early epoch. 



5. DISCUSSION AND CONCLUSION 

With our Chandra observation of SN 2002bo we have 
set the most sensitive X-ray flux upper limits at an ear- 
lier epoch than for any previous SN la. The previous 
best case was that of SN 1992A for which a ROSAT 
upper limit of ~10~^'* ergs cm~^ s""'^ was set on the 
0.5-2 keV band X-ray flux ~35 days after explosion 
(jSchlegel &: Petr ^119931 ). Our Chandra limit in the same 
energy band is an order of magnitude lower and was set 
only 9 days after explosion. We also set a sensitive upper 
limit in the important 2-10 keV band. 

Converting our flux upper limits to constraints on the 
density of the CSM in the system depends on uncertain 
assumptions about the thermodynamic state of the hot 
plasma in the expanding ejecta and shocked wind. As 
we have shown in this paper there is roughly an or- 
der of magnitude difference in the inferred wind den- 
sity under the assumption of fully equilibrated electron 
and ion temperatures vs. the case with non-equilibrated 
temperatures. The latter case yields the more conserva- 
tive (i.e., higher) constraint: a wind density parameter 
w < 1.2 X 10^^ g cm^^. In terms of a slow wind with 
velocity of Vwio — 10 km this corresponds to a up- 
per limit on the mass-loss rate of Af < 2 x 10~^ Mq 
yr~^. The much higher sensitivity of the Chandra data 
notwithstanding, this value is larger by a factor of 10 or 
so than the mass loss upper lim it derived from ROSAT 
data bv lSchlegel fc Petrd (fl993) for SN 1992 A. This dis- 
crepancy is due t o the simplicity of the model used by 
ISchlegei &: Petrd (|1993D to evaluate their upper limit: 
they do not include absorption by residual wind mate- 
rial above the forward shock, assume electron-ion tem- 
perature equilibration, compare to model luminosities 
without applying a bolometric correction, and calcu- 
late the SN age from the time of maximum rather than 
from the time of explosion. In other comparisons our 
Chandra upper limit is comparable to those found previ- 
ously from limits on the H ct flux for the n ormal type 
la supernovae SN 1994D (|Cumming et all [1996 ) and 
SN 2001el (jMattila et al.ll2005[ ). In general limits on M 
of '^lO"^ Mq yi:~^ are not stringent enough to rule out 
the class of symbiotic-type binaries as SN la progenitors 
at least not for these particular cases. 

The value of these results lies in our ability to calcu- 
late, using well understood physics, the expected X-ray 
emission from hot gas. Viewed in this light, we briefly 
discuss limits on the CSM obtained from radio nondetec - 
tions of a number of nearby SNe la ( Panagia et al.l l2006) . 
The radio results rely on semi-empirical parameterized 
functional forms for the time- and frequency-dependence 
of synchrotron emission and free-free absorption, whose 
relevant parameters are assumed to have values given by 
radio results for SNe Ib/c. Likewise the essential pa- 
rameter, i.e., the one linking the wind density parameter 
to the radio luminosity of the SN, cannot be calculated 
from theory with any accuracy, and an empirical calibra- 
tion, again drawn from measurements of SNe Ib/c, must 
be used. Although the sensitive radio flux limits clearly 
argue for low density environments around SN la, the 
extremely l ow numerical l imits on the mass-loss rates 
claimed by iPanagia et al.l (|2006( ) cannot yet be consid- 
ered definitive. 

The other normal SN la we study here is SN 2002ke. 
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We re-examine the claim bv llmmler et al.l (|2006f) of a ten- 
tative X-ray detection by Swift and find that we cannot 
substantiate it. We pay particular attention to the as- 
trometric and photometric calibration of the Swift X-ray 
data by comparing to a Chandra observation done several 
months later. We find no evidence for a significant X-ray 
detection of SN 2005ke by either Swift or Chandra, to a 
flux limit that is several factors higher than what we ob- 
tained for SN 2002bo. Since this limit is at a much later 
epoch, when the intensity of the CSM/ejecta interaction 
should be much reduced, we did not attempt to deter- 
mine numerical limits to the wind density parameter for 
this SN. 

We have also presented Chandra upper limits for the 
two examples of SNe la with clear evidence for circum- 
stellar interaction: SN 2002ic and SN 2005gj. The upper 
limit on the X-ray luminosity of SN 2002ic in the 0.5 — 6 
keV band unexpectedly reveals a serious drawback to 
the interaction model proposed previously (CCL): the 
predicted X-ray flux turns out to be larger by at least a 
factor of four. We identified the major missing element of 
the model responsible for the controversy: macroscopic 
mixing of cool metal-rich ejecta fragments into the hot 
gas of the forward shock, which results in strong absorp- 
tion of the X-rays emitted by the forward shock. In- 
terestingly, the absorption of X-rays by mixed shocked 
ejecta should have the effect of decreasing the required 
CS density in the model to explain the late time optical 
luminosity. This effect together with the higher radia- 
tion efficiency of the interaction with a clumpy CS mat- 
ter should result in a slightly higher expansion velocity 
of the shocked SN ejecta in the interaction model. These 
outcomes appear to be preferred by the optical spectra 
of SN 2002ic (see CCL). 

SN 2005gj appears to belong to the SN 2002ic-li ke fam- 
ily o f SN la with dens e CS envelopes ( Aldcring et al.l 
l2006HPrieto et al.ll2007t l. The Chandra upper limit from 
this object at about day 80 is a factor of four larger than 
the flux predicted by the interaction model of SN 2002ic 
recomputed for the corresponding epoch. Although the 
strength of the interaction argues against it, one possible 
explanation is that the CS density around SN 2005gj is 
just somewhat lower than around SN 2002ic. A better 
possibility invokes mixing of the shocked ejecta, which 
can reduce the emergent X-ray flux. Which of these is 
more likely to be the correct explanation requires a bet- 
ter understanding of the environment of SN 2005gj than 
we have at present. Future studies of the optical light 
curve and spectra of this interesting object are strongy 
encouraged. 

Resolving the issue of the X-ray non-detection of 



SN 2002ic-like objects has its dark side due to our in- 
troduction of the parameter Tqc, which is essentially in- 
calculable. Even three-dimensional hydrodynamical sim- 
ulations are unlikely to be able to recover this value in 
a fully self-consistent way. We, therefore, can predict 
in detail neither the spectrum nor the evolution of the 
X-ray flux from SN 2002ic-like supernovae. Two rele- 
vant remarks can be made, however. First, in the case 
of SN 2002ic the model without mixing predicts an in- 
creasing X-ray flux in the band below 10 keV up until day 
^400 (CCL). Therefore, X-ray detection at earlier epochs 
is unlikely to be more favorable, as the lack of detection 
of SN 2005gj at roughly 80 days after explosion (versus 
260 days for SN 2002ic) tends to support. The second 
remark relates to the spectrum of the emergent X-ray 
emission: hard X-rays (i.e., with photon energies greater 
than ~20 keV) are not affected by the absorption that 
mutes the lower energy flux (see Fig. [T2)) . Future sen- 
sitive X-ray observations covering a wider energy band 
(1 - 30 keV), therefore, should reveal a SN 2002ic-hke 
event as a strongly absorbed hard X-ray source and thus 
verify the proposed mixing scenario. 
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APPENDIX 

X-RAY EMISSION FROM HOST GALAXIES 

Only the host galaxies of the two nearby SNe, NGC 3170 and NCG 1371, show any evidence for X-ray emission. 
In each case there is a compact, spectrally hard, nuclear X-ray source in addition to faint diffuse emission. Here we 
present a brief report on the new information from the Chandra observations of these galaxies. 

NCC 3190 and other nearby galaxies 

NGC 3190 is a member of the Hickson Compact Group 44, also the Leo HI group, and is classifled as a Low Ion i zation 
Nuclear Emission Region (LINER) galaxy. X-ray emission from this galaxy was first reported by Pildis et al.' (19951) 
using a short (5 ks) i?OS'^TPSPC observation from which it was only possible to determine a broadband luminosity 
of 4 X 10^^ ergs s~^ (0.07-3 keV band), where we have corrected their published result for our assumed distance to the 
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galaxy. iLiu &: Bregmanl ()2005[ ) report the presence of a compact nuclear X-ray source with an estimated X-ray flux of 
8.8 X 10^9 ergs s"! (0.3-8 keV band) based on ROSATRRl data. 

The Chandra data (Fig. [T]) reveal that the compact nuclear source is spectrally hard and, in addition, is embedded 
in a faint, diffuse, spectrally soft, extended component. This latter component tends to lie north and east of the 
nucleus, and clearly avoids the prominent dust lane that runs through the galaxy from southeast to northwest some 
10" south of the nucleus. The spectra shown in Fig. [5] were extracted from the 5 regions shown on Fig. [TJ We have 
already discussed the bottommost spectrum in this figure (i.e., the background spectrum) in §3. All fits were carried 
out with the likelihood figure-of-merit function used for the background model fits, which does not yield an explicit 
goodness-of-fit criterion. However, the models, at least visually, provide a fairly good description of the spectral data. 
All spectral fits, except the absorbed nuclear power-law component, include absorption fixed at the Galactic HI column 
density value of — 2.1 x lO^*' cm~^. 

The top spectrum comes from the innermost region: a circle 1.5" (160 pc) in radius centered on the mean position 
of the hard nuclear source (R. A. =10:18:05. 63, Decl. =-1-21:49:56). For an acceptable fit the spectrum requires two 
components: a highly absorbed power-law as well as thermal emission (using the "mekal" model in xspec). The best 
fit thermal emission model has kT = 0.76 ±0.05 keV and emission measure, nennV = 4.3 x 10^^ cm^'^. The power-law 
component has a fixed photon index of F = 1.4 with a best-fit intrinsic absorption of A^h = 1-7 ± 0.3 x 10^"^ cm^^ and 
flux density at 1 keV of Fe = 6.2 ± 1.4 x 10"^ photons keV'^ cm'^ s^^. The unabsorbed flux (0.5-10 keV) of the 
power-law component is 5.5 x 10"^'^ ergs cm^^ s~^ or a luminosity of 3 x lO"'" ergs s~^. It is clear that the LINER 
nature of this galaxy is due to an obscured, low luminosity, active galactic nucleus (AGN). 

The three spectra in the middle of Fig. [5] are of the spectrally soft, extended emission and all are well described 
by thermal emission alone. The red spectrum was extracted from the second smallest region on Fig. [TJ an ellipse 
with semi-major axis lengths of 3.9" and 3.2" (420 pc x 340 pc). The best fit values for the thermal emission were 
kT = 0.80 ± 0.05 keV and nennV = 3.5 x 10^^ cm~^. The green spectrum came from the next larger region: an 
ellipse with semi-major axis lengths of 9.6" and 8.4" (1.0 kpc x 0.9 kpc). This was well fit by thermal emission alone 
with a best fit kT — 0.74 ± 0.08 keV and nennV = 2.9 x 10^^ cm~^. Finally the blue spectrum came from the next 
region: an ellipse with semi-major axis lengths of 45" and 25" (4.8 kpc x 2.7 kpc). Three obvious point sources were 
excluded from this region. This was well fit by thermal emission alone with a best fit kT = 0.37 ± 0.06 keV and 
emission measure, neUnV = 4.2 x lO^'^ cm~^. 

From the spectral extraction regions we estimate emitting volumes (assuming an ellipsoidal geometry with the line- 
of-sight depth equal to the mean of the axis lengths of the extraction regions) and then convert the emission measures 
given above into the density of the emitting plasma (for n^/nH — 1.2). Wc find density values that vary from nu — 0.27 
cm~^ near the nucleus through values of 0.069 cm~'^ and 0.015 cm~'^ for the inner and mid galaxy regions to a value 
of 0.0025 cm~^ for the outer galaxy. These values are broadly consistent with an r~^ profile perhaps pointing toward 
an outfiow from the central regions of the galaxy. In summary, th e X-ra y emission properties of NGC 3190 closely 
resemble those of other LINER galaxies studied bv lTerashima et al] (j2002[ ): a low luminosity and obscured AGN with 
soft [kT ^ 0.8 keV) extended thermal emission. 

X-ra y emission from two other galaxies (NGC 3193 and NCG 3185) in HCG 44 were also reported bv iPildis et al.l 
(|1995D . We also detect X -ray emission from these same galaxies in the Chandra data, which, thanks to a fortuitous 
value of the roll angle, happen to lie within the field of view. Although clearly detected, the emission from NGC 3193 
falls near the gap between chips S3 and S4, which renders derived results somewhat inaccurate. We do not consider 
this galaxy further. 

NCG 3185 is '--^9' off-axis where the imaging quality of Chandra is modest (^^5"). The PSF there is however sufficient 
to separate the comp act nuclear emission from the fainter off-nuclear source, XMMU J101737. 4-1-214144, discovered by 
iFoschini et al.l ()2002l ). The pubfished XMM-Newton luminosity for the off-nuclear source is 1.3 x 10^^ ergs s~^ (0.5-10 
keV band) assuming a power-law source with a photon index of 2.0 and distance of 21.3 Mpc. We detect this source in 
the Chandra data with a net total of 21 ± 5 counts above background from which we estimate a luminosity of 7 x 10^^ 
ergs s~^ (in the same spectral energy band using the same spectral form as the XMM-Newton study). This suggests 
some variability in the X-ray luminosity of this off-nuclear source. The nuclear source is stronger (detected at 72 ± 9 
counts above background) and generates a l uminosity of 2 xlOf^ ergs (again for the same band and spectral form). 
This agrees with the Lx value published bv lCappi et al.l ()2006l ) from XMM-Newton observations. 

NGC 1371 

NGC1371 has typically been classified (see, e.g., lEskridge et al.l [200^ as a weak bar, early-type spiral galaxy 
(SAB(rs)a according to the RC3). It is known to display an "E xtended Nuclear Emission- line Region" (ENER) 
in the form of a disk-shaped zone of Ha emission (jHameed &: Dever eux 1999) extending over an elliptical region with 
major axis lengths of roughly 30" x 13" (2.5 kpc x 1.1 kpc) (see Fig. |9]). A compact nuclear source appears in Ha 
as well. The total Ha luminosity of the galaxy is Lrq = 4.9 x 10""^ ergs s~^, which is uncorrected for extinction and 
therefore a luminosity lower limit. The nuclear component (i.e., the emission within 1 kpc « 12", which contains 
much of the ENER emission) contains roughly 6% of the total Ha luminosity. The nucleus of this galaxy has not been 
spectroscopically identified. 

Our Chandra data reveal a compact source of X-ray emission at R. A. =03:35:01. 35, Decl. =—24:55:59.6 that is coinci- 
dent with the nucleus and is well described by an absorbed power law spectrum (see top curve and data points in Fig. [6]). 
The best- fit parameter values are a photon index of F = 1.5^q'2, intrinsic absorption of A'h = 2.5^q'4 x lO^'^ cm~^ 
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(considerably higher than the Galactic column toward NGC 1371), and flux density at 1 keV of Fe — l.G^to's x 10 
photons keV^^ cm^-^ s^^. The unabsorbed flux (0.5-10 keV) is 1.2 x 10"^^ ergs cm~^ corresponding to a luminosity 
of 4 X 10^° ergs s~^. Whether the excess absorption is intrinsic to the nuclear source or comes from the intervening 
foreground portion of the galaxy is not clear. The nuclear X-ray and Ha luminosi ty values f all rig ht on the published 
correlation between these quantities for Seyfert galaxies, LINERs, and AGN from iHo et al.l (|2001l i. Thus we conclude 
that NGC 1371 contains a low luminosity AGN with properties similar to those found in other nearby spiral galaxies. 

We also detect diffuse X-ray emission on larger scales. The red spectrum in Fig. [6] (second from top) comes from the 
inner elliptical region identified on Fig. 31 which corresponds well to the extended Ha disk. The spectral extraction 
region used here has major axis lengths of 31.4" x 16.1" (2.6 kpc x 1.3 kpc) in size. Fits to the X-ray spectrum from 
this region prefer a two component model: low temperature thermal emission {kT « 0.3 keV) in addition to either 
a hot thermal [kT k, 10 keV) or power law (F « 1.8) component. Note that for these fits the absorbing column was 
fixed to the Galactic value (i.e., no extra absorption intrinsic to NGC 1371 was required). The unabsorbed 0.5-10 
keV band flux of the hard component alone (for either spectral form) is 4.0 x lO"^** ergs cm~^ s~^ corresponding to 
a luminosity of 1.4 x 10^^ ergs s~^. The soft component yields a flux of 1.6 x 10^^"* ergs cm^^ s^^ (equivalent to a 
luminosity of 5.6 x 10^^ ergs s""'^) entirely in the 0.5-2 keV band. 

Faint diffuse soft X-ray emission is visible on Fig. [3] within the outer elliptical region centered on the nucleus. This 
region, with major axis lengths of 104" x 75" (8.6 kpc x 6.2 kpc), covers the bright central portion of NGC 1371. 
Fits to the X-ray spectrum (shown as the blue curve and data points third from the top on Fig. [6|) are consistent 
with a thermal spectrum {kT « 0.4 keV). The X-ray flux is 2.8 x 10"^'' ergs cm~^ s""'^ (equivalent to a luminosity of 
1.0 X 10'^^ ergs s~^) and again appears entirely in the 0.5-2 keV band. 

The fitted emission measure of the soft thermal component in the disk region is nennV ^ 3 x 10^^ cm^'^, while 
that of the outer elliptical region is 4 x 10^^ cm~^. These values lead to mean gas densities of nn = 0.015 cm"'^ in 
the disk region and nn — 0.0024 cm~^ within the outer ellipse. For the inner disk- like region of NGC 1371, the close 
morphological match between the X-ray and Ha emission argues that the X-rays we see surrounding the nuclear AGN 
come from the integrated flux of sources linked to recent star formation, such as supernova remnants, winds from 
young massive stars. X-ray binaries, and so on. The composite nature of the X-ray spectrum, especially the presence 
of the hard emission component, points to the need for a contribution from a population of X-ray binaries. However no 
individual X-ray source is clearly visible in this region (aside from the nucleus) with a luminosity greater than ~10'^^ 
ergs s^"'^. 

There are three unresolved X- ray sourc e s that are nearly coincident (i.e., offset by <3") with individual H H regions 
in the Ha image of ,Hameed fc Devereuxl ()1999f ). The brightest of these, CXOU J033506. 0-245713, is approximately 
1.6' southeast of the nucleus and is indicated by the cross in Fig. [9] (right panel). The 0.5-10 keV band flux of this 
source is 3 x 10~^^ ergs cm~^ s"'^ assuming a hard power-law or high temperature bremsstrahlung spectral form, either 
of which describes the data weU. The other two sources, CXOU J033500. 2-245437 and CXOU J033452. 7-245527, are 
each roughly one-third as bright in the Chandra data. The luminosities of these sources, if they are in NCG 1371 as 
seems likely, are 3 x 10'^* ergs s^^ for the two fainter ones and 10"^^ ergs s^^ for the brighter one. With the exception 
of the nucleus, these are the brightest sources within the optical extent (5.6' x 3.9') of NGC 1371. Their high X-ray 
luminosity, generally hard spectra, and spatial association with H H regions argue that these are accreting high mass 
X-ray binaries. 

Four other sources were detected at high significance (>10 counts) within this area, while there were three other s with 
7-8 c ounts. None of these other sour ces shows a counter part in the Ha or continuum images of (Hameed & D evereuxl 
fl999l) or in 2MASS infrared images (jJarrett et al.l [20031 ). We expect only ~3 unrelated background sources at these 
flux levels, so most of these should be new sources associated with NCG 1371. 
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